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Abstract We explore an accretion model for low lumi- 
nosity AGN (LLAGN) that attributes the low radiative 
output to a low mass accretion rate, M a , rather than 
a low radiative efficiency. In this model, electrons are 
assumed to drain energy from the ions as a result of 
collisionless plasma microinstabilities. Consequently, the 
accreting gas collapses to form a geometrically thin disk 
at small radii and is able to cool before reaching the 
black hole. The accretion disk is not a standard disk, 
however, because the radial disk structure is modified by 
a magnetic torque which drives a jet and which is pri- 
marily responsible for angular momentum transport. We 
also include relativistic effects. We apply this model to 
the well known LLAGN M87 and calculate the combined 
disk-jet steady-state broadband spectrum. A comparison 
between predicted and observed spectra indicates that 
M87 may be a maximally spinning black hole accreting 
at a rate of ~ 10~ 3 M Q yr -1 . This is about 6 orders of 
magnitude below the Eddington rate for the same radia- 
tive efficiency. Furthermore, the total jet power inferred 
by our model is in remarkably good agreement with the 
value independently deduced from observations of the 
M87 jet on kiloparsec scales. 
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1 Introduction 

We summarise an accretion mod el for LLAGN that has 
previously been applied to Sgr A* ( Jolley fe Kuncid . l2007h . 
In this model, described in more detail below, initially 
collisionless gas accretes at a very low rate and collapses 
into a geometrically thin disk at small radii as a result of 
wave-particle resonances that facilitate efficient electron- 
ion coupling. Magnetic coupling between the disk and 
a jet is self-consistently modelled by a magnetic torque 
with a prescribed radial profile. Here, we apply this model 
to M87. 

The giant elliptical galaxy M87 (NGC 4486) is a 
LLAGN situated a t a distance of d = (16 ± 1.2) Mpc 
( Tonrv et all 1200 it ) in the Virgo cluster. It harbours a 
central supermassive black hole (SMBH) of mass M = 
(3.2±0.9) x 10 9 M Q , with a rapidly rotatin g disk of ionised 
gas, c onsisent with a keplerian thin disk (jMacchetto et all 
Il997l ). ac companying a prominent one-sided jet first de- 
tected bv lCurtisl (Il918l). The nucleus has a lum inosity of 



10 42 ergss" 1 (|Biretta. Stern fc Harris! . Il99lh which is 



at least two orders of magnitude below the luminosity 
expected for a standard thin accretion disk accreting at 
the Bondi rate Mb = O.lM ^yr" 1 , as determined from 
Chandra X-ray observations (|Di Matteo et al. , 2003). The 
observed luminosity from the nucleus is likely to be less 
than the total accretion power, however, because a sig- 
nificant proportion of the available accretion energy is 
used to power the large observed jet, with a total ki- 
netic power estimated to be as large as 2 x 10 43 ergs s _1 



(|Revnolds et allll996l ). 

The inflowing plasma in LLAGN is collisionless be- 
cause the dimensionless mass accretion rate is so low 
(rh = L/Lftdd *C 1) that there insufficient time for elec- 
trons and ions to come into equipartition via two-body 



processes before reaching the event hor izon ([Shapiro. Lightman fc Eardle^ 
ll976tllchimaruLll977t : lRees et al.Lll982h . If other coupling 
processes are unable to equilibrate the electrons and ions 
within the inflow timescale and if the ions are preferen- 
tially heated by viscous dissipation of the gravitational 
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binding energy, then the resulting accretion flow can- 
not radiate its internal energy before reaching the hole. 
This lea ds to a Radiatively Inefficient Accretion Flow 
(RIAF) (jNaravan fc Yil . ll994l ). RIAF models for LLAGN 
attribute the low luminosity to a low radiative efficiency 
e, rather than a low mass accretion rate, M a . Recent 
variations on the basic RIAF model consider a reduced 
accretion rate close to the black hole du e to convec- 
tive motions ( Quataert fc Gruzinovl . I2000T ) or outflows 
( Blandford fc Begelmanl I1999D . A reduced M a at small 
radii appears to be necessary, at least for Sgr A* , where 
polarization mea s urem ents imply M a < 4 x 10 _8 M Q yr^ 1 
(|Macauart et al.[|2006l) . 

M87 has a large (« 2 kpc long) relativistic jet. Unlike 
winds, relativistic jets cannot carry away large amounts 
of mass, and therefore cannot be responsible for reducing 
the accretion rate close to the black hole. Furthermore, 
the mechanism by which the jet is coupled to the un- 
derlying accretion flow has yet to be explicitly modelled, 
so the effect of the jet on the radial structure of the 
underlying accretion flow is not known. Here, we con- 
sider the possibility that the inflowing gas in M87 and 
other LLAGN is accreting at a very low rate and cools 
and settles into a compact, geometrically thin disk at 
small radii. We suggest that the collisionless accreting 
plasma can form a cool disk because plasma microinsta- 
bilitics act immediately to reduce any large temperature 
gradient between the electrons and ions. Wave-particle 
resonances can the n efficiently couple the elect rons and 
ions, as proposed by Begelman fc Chiuehl (fl988l) fsee also 



A discussion of this work and some concluding remarks 
are given in Section 



2 Conditions for a Cool Disk 

Here, we describe the properties of collisionless accretion 
flows when the assumption of a two-temperature plasma 
ceases to remain valid. Th e deta iled derivations are pre- 
sented in lJollev fc Kuncid (|2007l ). 

If there exists a mechanism to transfer internal en- 
ergy from the ions to the electrons on an inflow time 
(iBegelman fc Chiuetj.ll988tlBisnovatyi-Kogan fc Lovelace . 
| 200(ilQuataerd.ll998l:lGruzinovl . ll998l : lQuataert fc Gruzinovl . 
Il999l; iBlackmanl Il999f ).~ then the accretion flow geom- 
etry will deflate from an ion-pressure-supported, two- 
temperature torus predicted by RIAFs to a much thin- 
ner, quasi-thermal structure. For a thermal pressure sup- 
ported flow at the electron virial temperature the height- 
to-radius ratio is h/r = (rric/nip) 1 / 2 ps 0.02. This defines 
a geometrically thin disk. 

The condition that the ratio of the cooling timescale 
to the inflow timescale i c ° o1 /£ inflow < 1 is required for the 
disk to cool before t he gas reaches th e black hole. This 
implies (see also e.g. iRees et al 
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where e = O.leo.i is the radiative efficiency and 



Bisnovatvi-K ogan fc Lovelacd2 000: Quataert 1998; G rlMl 1 ""- 1 is the dimensionless v iscosity parameter fr omjgati- 
19981 : buataert fc Gru Z inoJl999riBlackmaJl99l . The dkrd accretl0n disk theor y (|Shakura fc Sunyaevj , 11973J). 



presence of cold molecular gas inside the Bondi radius 
( Tan et al 1 l2007h may indicate the presence of a cool, 
thin accretion disk. This is inconsistent with RIAF mod- 
els whic h form a geometr ically thick, hot ion torus ge- 
ometry ( Rees et all Il982l ). 

In what follows, we wish to determine whether the 
low luminosity of M87 (and other LLAGN) can be at- 
tributed to a low M a rather than a low radiative effi- 
ciency. We proceed by considering a geometrically thin, 
cool, single temperature relativistic accretion disk that is 
modified by magnetohydrodynamic (MHD) stresses. The 
disk is coupled to a relativistic jet via an MHD torque 
acting a cross the disk sur f ace. T his model is described in 
detail in lJollev fc Kuncicl ( 20071 ). The radiative efficiency 
of our cool disk is somewhat lower than that of the stan- 
dard Shakura-Sunyaev disk (jShakura fc Sunvaevl . Il973t ) 
as a result of efficient extraction of accretion power by 
the jet. In Section[2j we examine the physical conditions 
needed for a cool disk. In Section [31 we present a sum- 
mary of the relevant equations for the modified disk flux 
and the steady-state spectrum resulting from a jet mag- 
netically coupled to the underlying accretion flow. We 
compare the spectrum predicted by our coupled disk-jet 
model with the observed spectrum for M87 in Section [H 



Hence by relaxing the assumption that electrons and ions 
can interact only through Coulomb collisions, the colli- 
sionless accretion flow in low-m systems must be geomet- 
rically thin. The condition (TT]) implies that the collision- 
less accretion flow is able to cool on an inflow timescale 
and thus, is not advective. 



3 Coupling a Magnetized Jet to a Relativistic, 
Cool Disk 

3.1 Non-Zero Torque on the Disk Surface 

The MHD stresses present in an accretion disk can pro- 
duce a torque acting across the disk surface as well as 
at the last marginally stable orbit. The surface torque 
can efficiently remove angular momentum vertically out- 
wards from the disk and does work against the d isk, di- 

recting energy vertically to form a magnetized jet ( Kuncic fc Bicknelil . 
12004b . 

The radiative flux fr om a relativistic, torqued disk is 
dJollev fc Kuncid . l2007h 
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where x = r/r g , / NT (x) is the lNovikov fc Thornd (jl973h 
relativistic correction factor, f^(x) is a correction fac- 
tor f or a nonzero to r que (N ZT) at the inner disk bound- 
ary (| Agol fc KrolikL l2000h . and f^(x) is an analogous 
correction factor for a nonzero torque on the disk surface. 

Using global energy conse rvation the disk flux p ro- 
file ((2]) can be expressed (see IJollev fc Kuncid 120071 for 
details): 
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is the efficiency of the torque acting at the last marginally 
stable orbit, and 
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is the efficiency of the jet. The functi ons C(x) and A(x) 
are rel ativistic correction factors from lNovikov fc Thorn 



1973), and the following are derived in lJollev fc Kunci 
2007h : 
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The mass accretion rate can be written as 

M a = ^ (9) 

where Lgdd = 47rGMm p c/er T . 

The input parameters for the modified disk model 
are the dimensionless black hole spin parameter a, the 
fractional efficiency of the torque at the last marginally 
stable orbit Z\e/e NT — 0.10, and the fraction of accretion 
power injected into the jet, ej/e a - This last parameter 
has an upper limit in order for the disk flux to remain 
positive at all radii. 

The effect of the nonzero torque across the disk sur- 
face is to do work against the disk, thus reducing the disk 
flux over the range of radii where the magnetic torque 
is strongest. This counteracts the effect of the nonzero 
torque at the inner disk boundary, which enchances the 
disk flux The combined effects of these two 

torques is clearly evident in the radial flux profiles in 
Fig. [T] The resulting disk flux radial profiles are sub- 
stantially modified from their corresponding zero-torque 
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Fig. 1 Radial flux profiles for the jet-modified disk model 
with different model parameters: a is the black hole spin pa- 
rameter, e a is the overall accretion efficiency and ej/e a is the 
fractional jet power. The solid line corresponds to a relativis- 
tic disk torqued at the inner boundary and on its surface; 
the dotted line corresponds to the same disk without torque 
effects. 



profiles (Fig. [TJ dotted lines). It is clear that the nonzero 
magnetic torque acting on the disk surface results in a 
disk radiative efficiency £d that is lower than that of a 
non-torqued disk. 

We have explicitly taken into account how the local 
disk radial structure is modified by a magnetized jet that 
is primarily responsible for angular momentum trans- 
port. This results in a disk spectrum that is modified 
with respect to that predicted by standard theory. 



3.2 Jet Emission 

Wc identify the nonzero magnetic torque across the disk 
surface with the mechanism responsible for extracting 
accretion power from the disk and injecting it into a jet. 
Some of the magnetic energy is subsequently converted 
into kinetic energy. We expect a fraction of the parti- 
cles to be accelerated to nonthermal, relativistic ener- 
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gies. Synchrotron radiation by relativistic electrons will 
then contribute significantly to the radio emi s sion. 

Following the method in I.Tollev fc Kunticl (j2007h . we 
divide the jet into a series of quasi-cylindrical sections of 
thickness Az, and calculate the total emission spectrum 
by summing the contributions from each component. 

We consider a relativistic jet with bulk Lorentz factor 

J] and Doppler factor £ = {/] 1 - (1 - rr 2 ) 1 / 2 cos 0j j 

where 6j is the angle between our line of sight and the 
M87 jet axis. We use the following simple radiative trans- 
fer model to calculate the observed specific luminosity 
due to the net contribution from each jet component 
(assuming isotropic emission in the source rest frame): 
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where AA « TrrZizsin^j is the projected surface area of 
each emitting cy linder, £, s , y . n is the synchro tron source 
function (see e.g. iRvbicki fc Lightmanlll979l for relevant 
formulas) and 



r syn 



<rV, yn As 



is the synchrotron optical depth along a path length As 
through each cylindrical section. 

The electron number density N c and hence the mag- 
netic field B decline with jet height z according to 
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where the first term in square brackets refers to the 
bulk jet kinetic energy and the second term refers to 
the electron kinetic energy and the magnetic energy. The 
equipartition factor f cq = 1 is used to relate the mag- 
netic and electron energy densities. Equation (11) is used 
to calculate N c at the base of the jet. 



4 Application to M87 



Figure [2] shows the predicted combined disk and jet spec- 
tra for our model for different parameters. The obser- 
vational data p oints (corrected for extinction) are taken =? 
fromlHol (jl999l ) (plus signs) and lMeisenheimer. Roser fc Sc hlol 1 <^'bu 1 ri rame ers 
( 199a ) (diamonds). The inclination of the nuclear disk a £ j/ e 



Log [v (Hz)] 



Fig. 2 Observed and predict ed spectra for M87. The plus 
symbols are data points from Ho (1999) and the diamonds 
are from Meisenheimer et. al. (2007). The solid line is the 
theoretical disk+jet steady-state spectrum predicted by our 
model. The dashed line is the disk spectrum, and the dotted 
line is the synchrotron jet spectrum, (a) is for a black hole 
with zero spin, and (b) is for a maximally spinning black hole. 
The launching height and the bulk Lorentz factor for the jet 
are given by zo and Ij, respectively. See Table 1 for other 
parameters. 

Table 1 Parameters used in the disk model. The black hole 
spin is a, £j/e a is the maximum allowable fraction of the ac- 
cretion power removed by the jet, e a is the total accretion 
efficiency, ed is the disk radiative efficiency, M a is the accre- 
tion rate in M©yr -1 and Pj is the total jet power in erg s _1 . 
The black hole mass is M = 3.2 x 1O 9 M , the disk luminosity 
as a fraction of the Eddington luminosity is rh = 1.8 x 10~ 6 
and the optically-thin synchrotron spectral index is a — 0.99. 



inferred parameters 



rotat ion axis to our line of sight is 6& ss 60° (jMacchetto et akl g qq 
Il997l ). and the jet incli nation angle is 0\ ~ 30°, wit h an q.99 
ope ning ang le of « 60° (|Lv. Walker fc Junorl l2007h (see — — 
e.g. lBirettalll999l for a review of the properties of the M87 
jet). The disk luminosity as a fraction of the Eddington 
luminosity is id/iEdd = th = 1.8 x 10~ 6 . Table 1 lists 
the other physical parameters used in our model for the 
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zero spin (a = 0) and maximally spinning (a « 1) black 
hole cases. 

An accretion disk around a spinning black hole can 
extract more accretion power, P a = e a M a c 2 , than a disk 
around a non-spinning hole for the same accretion rate. 
This is because the last marginally stable orbit for a = 
0.99 is much smaller than that for a = (where r ms pa 
6r g ). Thus, the accretion disk in the high-spin case (with 
r ms pa 1.45r g for a — 0.99) reaches higher temperatures 
and emits a bluer spectrum than that in the zero spin 
case. 

Similarly, the jet launching height zq determines the 
normalization of the relativistic electron number den- 
sity N e and hence the jet synchrotron spectrum. Both zq 
and -Tj are constrained by the requirement that the jet 
remains non-dissipative (i.e. radiatively inefficient) and 
relativistic. These jet constraints produce a better overall 
agreement with the data for the high-spin case. Further- 
more, the jet power predicted for the maximally spinning 
(a = 0.99) case is in excellent ag reement with th e value 



There is observational evidence for the presence of 
cold gas in the vicinity of the nuclear disk in M87, in the 
form of molecular gas inside the Bondi radius. There has 
been speculation that the mass accretion rate could be 
reduced due to star formation in these regions, although 
any de finitive ev i dence for this scenario has yet to be 
found (|Tan et alj . [2001 . 



Pj ps 2 x lO^ergs" 1 deduced by iRevnolds et alj (|l996ft 
from observations of the M87 jet on kiloparsec scales. 
Our model thus predicts that the black hole in M87 may 
be maximally spinning. 



5 Discussion and Conclusion 

We have presented a model which combines existing the- 
ory for relativistic disk accretion with an explicit pre- 
scription for disk-jet coupling via a magnetic torque on 
the disk surface. The torque efficiently extracts angu- 
lar momentum and energy from the disk at small radii 
to drive a magnetized jet. Using this model, we have 
demonstrated that the low radiative output from M87 
can be attributed to a low mass accretion rate rather 
than a low radiative efficiency. From the predicted com- 
bined jet and disk spectra, our model indicates that M87 
may be a rapidly spinning black hole with a dimen- 
sionless spin a pa 0.99. We predict a mass accretion 
rate M a pa 1 X 10~ 3 M Q yr -1 and a disk radiative effi- 
ciency ed pa 0.01. This interpretation of the nature of 
black hole accretion in M87 differs from that of radia- 
tively inefficient models, which attribute the low lumi- 
nosity to an unusually low radiative efficiency, typically 
< 10~ 5 . Whereas radiatively inefficent accretion flows as- 
sume a thick, bloated torus geometry with a hot two- 
temperature plasma, our model assumes a geometrically 
thin, quasi-thermal disk. 

Our model for low-m accretors requires an efficient 
thermal coupling mechanism between the electrons and 
ions facilitated by collisionless plasma instabilities, re- 
sulting in a geometrically thin disk which can cool on 
an inflow timescale. The exact nature of the collisionless 
plasma instabilities required is currently the subject of 
future work. 
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